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Abstract 

The possibility of the determination of the neutron mean square charge radius from high-precision 
thermal-neutron measurements of the nuclear scattering length and of the scattering amplitudes 
of Bragg reflections is considered. Making use of the same interferometric technique as Shull in 
1968, the scattering amplitudes of about eight higher-order Bragg reflections in silicon could be 
measured without contamination problem. This would provide a value of the neutron charge radius 
as precise as the disagreeing Argonne-Garching and Dubna values, as well as a Debye- Waller factor 
of silicon ten times more precise than presently available. 
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The so-called mean square charge radius of the neutron (r^) is an important structure 
constant which reflects the internal charge structure of the neutron. In terms of the generally 
used form factors, (r^) is related to the derivative of the Sachs form factor G% at vanishing 
transferred momentum Q, 

1 dG^ 



(1) 
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Low-energy neutron-atom scattering experiments have been found useful |1|, |2[ to determine 
(r^) via high-precision measurements of the neutron-electron scattering length 6ne; 

^•4^(0. (2) 

where is the mass of the neutron. Because of the fundamental importance of these 
quantities for our understanding of the nucleon many ambitious experiments have been 
performed in the last decades (see e.g. and references therein). 

Despite these efforts the determination of 6ne is still unsatisfactory because there exist 
two sets of results which differ more than three standard deviations from each other. In this 
paper we propose an independent method based on high-precision measurements of neutron 
Bragg reflections on silicon. Assuming for all higher-order reflections a precision similar to 
that achieved by Shull [§] one can solve the discrepancy in the values of 6ne- 

In a flrst theoretical estimate, 6ne is expected to be given by the Foldy scattering length 

i, 

^aieory _ _i.46797i(4) X 10-^fm, (3) 

corresponding to the well-known value of the anomalous magnetic moment of the neutron 
0. The actual value can be deduced from low-energy neutron-atom scattering, for which 
the "scattering length" reads [|l|] 

h{Q) = 6nuclcar-&nc^[l-/(Q)] (4) 

" V ' 

electrostatic 

= ^nucleus + KcZf{Q) , (5) 

electrons 

where fcnuciear is the nuclear-interaction scattering length, ^nucleus is the neutron-nucleus scat- 
tering length (nuclear and electrostatic interactions), Z is the atomic number and f{Q) 
is the atomic form factor normalized to the forward direction: /(O) = 1. In Eq. (^), the 



Q-dependent electrostatic term is typically three orders of magnitude smaller than the con- 
stant nuclear term, hence it is reasonable to speak of "scattering length" despite the Q 
dependence. To get a sufficiently high sensitivity to this term scattering processes at Q 7^ 
(hence < f{Q) < 1) and high Z should be chosen. 

The best experimental values of b^^ (see e.g. Table I in |]^) can be grouped into two 
sets. Those deduced from thermal-neutron angular scattering on noble gases and from 
thermal- and epithermal-neutron transmission on liquid lead and bismuth [§, ^ agree on the 
value = —1.31(3) x lO^'^ fm (Argonne-G arching). The second set of 6nc values obtained 
from another transmission experiment on Bi |lO| and a measurement of thermal-neutron 
Bragg reflections on tungsten monocrystals [|Tl| are centered about b^^ = —1.59(4) x 10^^ 



fm (Dubna). There is a clear discrepancy because the 6ne values of the two sets differ more 
than three standard deviations from each other. In addition they differ about 10% from 
^ne^°^^) one being smaller than b'^^"^^, the other larger. While the difference in the 6ne values 
obtained from the Bi-transmission experiments |l^] is well understood ||^, ^ |12|, the 
Bragg-reflection method of Ref. ||lT| has never been revisited. In the latter experiment, the 



value of 6ne is extracted from the slope of b{Q) as a function of 1 — f{Q) [see Eq. (^]. The 
fitted data are the forward value fenuciear? obtained from a Christiansen-filter experiment, and 
the values for eight Bragg refiections b{Qhki)- The transferred momentum Qhki reads 



Qhki/^T^ = VWT¥TP/2ao = smOhki/K (6) 

where (hkl) are the Miller indices, is the side of the conventional cubic unit cell, O^ki is 
the Bragg angle and A is the neutron wavelength. 

For two different crystals, two very different values of 6ne were obtained in Ref. [|ll| : 



i.e., bne = —1.06 and —2.2 x 10^^ fm, which led the authors to postulate the existence of 
an additional scattering process. We rather consider that this discrepancy reveals the sys- 
tematic uncertainty of the experiment, where b{Qhki) is deduced from an integral-intensity 
measurement of the Bragg peak. This systematic uncertainty has probably been underes- 
timated in Ref. [|11[| for different reasons: (i) extinction is neglected, (ii) the uncertainty 



on the temperature factor (which is shown below to be of crucial importance) is neglected. 
Consequently, the accuracies obtained for b{Qhki) (typically 0.01 fm) seem unrealistic as 
compared with accuracies generally obtained by this kind of measurement (typically 0.1 fm) 
|T^. A hint to our conjecture is also given by the strong disagreement between other results 
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Tj] obtained inter alia with an integral-intensity measurement on monocrystals, and the 
Christiansen-filter results of Ref. [jT5[ for W isotopes. 

Let us now consider the method proposed by Shull in Ref. [Q, which also aims at mea- 
suring b{Qhj:i) for Bragg reflections on monocrystals. We refer to Ref. for details about 
the experimental setup. A collimated full-spectrum incident beam is Bragg reflected on a 
monocrystal blade (typical thickness t = 1 cm) in Laue transmission geometry. The wave- 
length selection is made through the Bragg condition (P), the crystal and the detector being 
moved simultaneously to form angles 6 and 26 respectively with the incident-beam direction. 
As shown in the dynamical theory of Bragg reflection for neutrons P, |l^ , the reflected beam 
inside the crystal consists of two coherent waves, the so-called Pendelldsung (pendulum so- 
lutions), interfering with each other and generating Pendellosung oscillations. An entrance 
and a scanning gadolinium slits (width 0.13 mm) are placed on both faces of the crystal to 
extract these fringes. Making use of the symmetry of the interference pattern a very precise 
alignment of the scanning slit with respect to the entrance slit is possible. The intensity at 
the center of the interference pattern is then measured as a function of the wavelength and 
compared with the theoretical expression. This expression may be obtained either by the 
Kato spherical- wave dynamical theory |16[ or by solving the Takagi-Taupin equations [H. 



For a perfectly-flat crystal and with infinitely-narrow slits, this intensity reads 

hki{\) oc h{\)AF^u?Jl { 3 ^'^r'.u ^l ' (7) 

where /q is the incident intensity, Jq is a Bessel function and \Fhki\ is the unit-cell structure 
factor of the crystal. This expression has to be slightly modified to take into account the 
finite curvature of the crystal, which adds the curvature radius R (typically 20 km) as a 
parameter to be fitted to the data, and to take into account the finite opening of the slits 
||T6| . Typically, forty fringes are observed for neutrons ranging the full spectrum of a thermal 
reactor. The structure factor can then be deduced with high precision from the period of 
these fringes. ShuU's method is thus in essence an interferometric method, which explains 
why it provides very high accuracies (contrary to the integral-intensity method). 

For a diamond-structure crystal like silicon, the unit cell consists of four elementary cells 
with two atoms. The corresponding structure factor reads [O 



hkl 



4X + X 6^eas (Q/.fc/) , (8) 
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where 6meas(Q) is related to b{Q) through the Debye- Waller temperature factor |T3 



bmeUQ) = KQ) X exp [-B{Q/Any] . (9) 

To get a precise b{Qhki) value, as needed to estimate fene, one needs precise values of both the 
structure factor F^ki and the temperature parameter B. For instance, for the (111) reflection 
on Si {Z = 14, ao = 5.43072 A), Shull's method provides pf b^eUQiu) = 4.1053(8) fm. 



Using the Debye- Waller factor obtained by measuring X-ray Pendellosung fringes at room 
temperature [|17| 



E = 0.4613(27) (10) 

one gets 

6(giii) = 4.1538(11) fm, (11) 

where three additional units in the error are due to the error on B. Equation (^ shows that 
this additional uncertainty increases exponentially with (see Fig. |I|). 

Following the approach of Ref. ||TB|, we can deduce ^nuclear from value (0) with Eq. 
(I), using the Argonne-G arching b^c and the atomic form factor DTSf: /(Qui) = 0.7526 (a 
precise value of f{Q) is actually not required 0]). This yields &nuciear = 4.1495(11) fm (which 
slightly differs from the value of Ref. |16| because of the more recent temperature factor), in 



agreement with the more precise value obtained recently by non-dispersive-sweep neutron 
interferometry 



&nuclear= 4.1507(2) fm. (12) 

Conversely, b^e can be deduced from values (|Tl|) and (|1^) with Eq. (|^), which provides (see 
Fig.0) 

bf^ = -0.89(32) X 10^3 fm. (13) 

This value agrees with the Argonne-Garching value but does not exclude the Dubna value 
(see Fig. ^). The same method applied to germanium (diamond structure with = 5.6575 
A) using = 8.1929(17) fm |2|, 6^eas(Qiii) = 8.0829(15) fm [0], = 0.57(1) A^ 



H and /(Qui) = 0.8542 |T|, leads to the value b^^ = 0.28(83) x 10"^ fm. Ahhough 



measurements on Ge (Z=32) could in principle lead to a value of b^e two times more precise 
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than on Si {Z = 14), this advantage is compensated by the two-times-larger scattering 
length, which implies a two-times-larger absolute error on h{Qhki)i the relative accuracy of 
Shull's method being 0.02% in both cases. Moreover, no ^nuclear value for Ge of comparable 
precision to (|T^) is available at present. Hence we only consider Si in the following. 

In Fig. |I], values (|r^) and (JT^) are represented as a function of 1 — /(Q). It is seen that the 
(111) reflection is not the optimal choice to calculate the slope of the line since it is close to 
the origin. Higher-order reflections are more appropriate, as shown by the simulated data. 
The error bars of these simulated points are calculated with the error as = 0.0027 on 
the temperature factor ([To|) only, assuming an ideal experiment with an infinite accuracy 
on 6meas- With such large uncertainties, there is no hope to reach a high precision on b^e, 
in particular to distinguish between the Argonne-G arching (solid line in Fig. |lD and Dubna 
(broken line) values. However, if enough (at least two) reflections are measured, both B 
and 6ne could be deduced from the data. The precision of the extracted B and bne can be 
improved by increasing the number of reflections. 

Let us now determine the possibly-observable reflections. A diamond-structure crystal 
has a face-centered cubic lattice, which implies that all Miller indices have to be either even 
or odd [|13|. Moreover, for a diamond structure, the structure factor (|]) leads to three types 
of reflections: (i) forbidden when h + k + I = 2 + An, where n is natural, (ii) weak when 
h + k + l is odd, (iii) strong when h + k + I = A + An. A thermal-neutron beam has typically 
a maximum flux for A = 1.2 A and the accessible angular range of the setup is assumed to 
be < 2^^ < 110°. These conditions combined with Eq. (|^) imply that sixteen reflections 
between (111) and (642) could be measured with reasonable intensities for Si. However, if 
the incident beam has a typical full spectrum 0.8 < A < 2.5 A, contamination has to be 
taken into account. For instance, while the (111) reflection is pure for 15 < 26 < 45° thanks 
to the absence of the (222) reflection which is forbidden (hence the result of Refs. P, p!6[|), 
(333) and (444) would be measured on A5 < 26 < 110° and 61 < 2^ < 110° respectively. 
This mixing implies that the reflected intensity would not be described by Eq. (|^ and that 
no structure factor could easily be extracted. Among the sixteen considered reflections, 
seven could be contaminated. The remaining nine reflections are listed in Table |. Three 
crystal blades would suffice to measure them all [for instance, a crystal blade cut parallel to 
the (220) planes could be used to measure the (111), (422), (511), (533), (711) and (551) 
reflections] . Among these reflections, three are strong (underlined in Table |) ; the intensities 
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range between 0.7 and 1.7 times the intensity of the already- measured (111) weak reflection. 
Finally, Eq. (|^ allows to estimate the number of fringes which would be measured for each 
of them. For instance, 44 fringes would be observed for the (711) reflection if t = 1 cm. 

These considerations show that these higher reflections should not be more complicated 
to measure than the (111) reflection already measured in Refs. |rB|. Hence, we can 



reasonably assume that the precision on 6mcas would be equal to that of Ref. [|TB|: 0.0008 A. 
With this assumption, the attainable precision for B and 6ne can be estimated. For B, we 
apply the approximate linear relation [Eqs. (§) and @] 

ln6meas(Q) ~ In ^nuclear " B{Q/47lf, (14) 

while for 6^0 we use Eq. (p. The uncertainty of a linear-fit slope can be calculated from the 
abscissas of the experimental points and from their errors Measuring the three strong 
reflections would give 



ae = 0.00040 a^^^ = 0.11 x 10"^ fm, (15) 

while with the eight new reflections of Table |, one gets 

as = 0.00027 A^ = 0.06 x 10"^ fm. (16) 

The temperature factor would be ten times more precise than the most precise value (|10[) 
and the 6ne value would be nearly as precise as the Argonne-Garching and Dubna values 
(see Fig. H). Since the present method is similar to the Dubna method on W, such a new 
result would certainly be of great help to decide the b^e discrepancy. 

In conclusion, the measurement of Pendellosung fringes in Bragg reflections is a very 



powerful technique. This is known for X-rays [|l^, |2J] but generally ignored for neutrons. 
By confirming ShuU's pioneer conclusion 

". . . the high sensitivity of the fringe positions . . . can be exploited . . . [for] ex- 
ample ... in assessing the neutron-electron interaction strength. . . " 

the present work opens new prospects for the field of Neutron Pendellosung Interferometry. 

We acknowledge useful discussions with H. Ranch, H. Kaiser and E. Jericha. J.-M. S. 
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neutron optics and by the National Fund for Scientific Research (Belgium). He thanks the 
members of the Atominstitut for their hospitality. 
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FIG. 1: Experimental ^nuclear of Eq. ( [l^ ) and ^(Qiii) of Eq. ( [LlD compared with the theoretical 
curves (Q) for ^nuclear = 4.1507 fm, 6ne = —1.31 (Argonne-Garching) and —1.59 x 10~^ fm (Dubna). 
The simulated points correspond to the reflections of Table ||; they are calculated with the Argonne- 
Garching value, with error bars only due to the uncertainty on the temperature factor (|lO|). 
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FIG. 2: Comparison of the theoretical value @ with the Argonne-Garching and Dubna 
experimental values and with the result (|l^). The simulated point has the Argonne-Garching 
value with the error bar of Eq. (|l^). 
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TABLE I: Weak and strong (underlined) Bragg reflections which could be measured with a thermal- 
neutron reactor without contamination problem, (hkl) are the Miller indices, f{Qhki) is the atomic 



form factor |18], A the neutron-wavelength interval, 26 the corresponding angles between the in- 
cident and reflected beam 
intensity through Eq. (S). 



cident and reflected beams. is the squared structure factor which provides the reflection 



(hkl) 


/ 1 J. 1. 7 1 

J \^ riKt J 


A (A) 


2e(°) 


|F/,hP (A2) 


(111) 


0.7526 


0.8-2.5 


15-47 


540 


(422) 


0.4788 


0.8-1.8 


42-110 


918 


(511) 


0.4600 


0.8-1.7 


45-110 


448 


(531) 


0.4150 


0.8-1.5 


52-110 


421 


(620) 


0.3902 


0.8-1.4 


56-110 


811 


(533) 


0.3764 


0.8-1.4 


58-110 


396 


(551) 


0.3432 


0.8-1.2 


63-110 


372 


(711) 


0.3432 


0.8-1.2 


63-110 


372 


(642) 


0.3249 


0.8-1.2 


67-112 


715 
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